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Abstract. A method is proposed for the characterisation of the entropy of cellular structures, based on 
the compactivity concept for granular packings. Hamiltonian-like volume functions are constructed both 
in two and in three dimensions, enabling the identification of a phase space and making it possible to take 
account of geometrical correlations systematically. Case studies are presented for which explicit calculations 
of the mean vertex density and porosity fluctuations are given as functions of compactivity. The formalism 
applies equally well to two- and three-dimensional granular assemblies. 

PACS. 82.70.Rr Foams - 81.05.Rm Granular materials - 05.20.Gg Statistical physics 

1 Introduction mained underdeveloped for several reasons: (a) it has been 

unclear why equlibrium-based description should apply to 



^ Many foams and cellular systems appear to evolve into these non-equilibrium systems; (b) it has been difficult to 

C3 ; steady states of snriilar structures irrespective of the spe- identify a phase space that makes it possible to either dis- 

cific dynamics and of the lengthscales over which the dy- entangle the geometrical correlations or at least enable to 

namics take place. The lengthscales can differ by orders treat them systematically; (c) it has been unclear what 

of magnitude between systems p[|. This prompted the use should play the role of a Hamiltonian, which is at the 

of statistical mechanical methods to analyse some of the foundation of the conventional statistical mechanical for- 

characteristics of such structures 2 . This approach has re- mahsm; (d) the extension of the concept of temperature 



RB is grateful to Prof N. Rivier for an illuminating discus- is far from obvious in these athermal systems. 
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This paper proposes a formalism based on the concept grains that an assembly of many grains can take. To move 

of compactivity which resolves all these difficulties. In from configuration to configuration these athermal sys- 

the first part we adapt a basic formalism, originally pro- tems require external agitation. Such agitation has been 

posed for granular systems, to planar cellular systems 0: . quantified by a new parameter, compactivity, defined as 

A convenient volume function is proposed as the natural X = dV/dS. 

analog of the Hamiltonian in conventional thermodynamic In cellular systems the situation apears to be differ- 

systems. The volume function makes it possible to iden- ent. Not only the dynamics by which configurations may 

tify a phase space and a method is outlined to isolate and change are different but also foams, for example, appear 

systematically analyse the correlations due to topological to be always in non-equilibrium, with bubbles ever grow- 

constraints. In the second part we use the insight gained ing in size. Nevertheless, it is often observed |2 that foams 

from planar systems to construct the formalism for three- and cellular systems converge to a scaling regime that can 

dimensional structures. The key step is the identification be regarded as a steady state. In foams the distributions 

of an exact volume function which makes it possible to of several structural properties (e.g. the volumes of cells) 

pinpoint the configurational phase space and evaluate its is stationary when the space dimension is scaled by an ap- 

dimensionality. We illustrate the formalism with cxphcit propriate power of time, r r' — r/t^. The existence of 

calculations of several case studies. a steady state suggests that these out-of-equilibrium sys- 



The compactivity concept [5] has been introduced in 
the context of granular packings to enable the use of statis- 
tical mechanical methods for the analysis of vibrated gran- 
ular systems. Relying on observations that under speci- 
fied vibrations granular assembliess settle into packings 
of reproducible densities, it was proposed that there is a 
steady state distribution of configurations. By identifying 
the configurations a statistical mechanical-like approach where are the degrees of freedom that comprise the 



tems may be analysed using statistical mechanical meth- 
ods Following the same logic as in granular systems 
one can apply the concept of compactivity to foams and 
cellular systems. In this formalism the partition function 
of a canonical ensemble in d dimensions is 



e-^T^e({q„})J]dg„ (1) 



can be constructed. However, thermal energy, which is phase space and 0{{qn}) is the probability density of find- 
relevant in conventional molecular systems, is irrelevant ing the system at a particular point in this space. This 
in granular systems. This led to the suggestion to replace density is nontrivial due to the connectivity forced by the 
the conventional Hamiltonian by a volume function W. mechanical equilibrium. A is the analog of Boltzmann's 
The entropy of such systems is volumteric and it depends constant and /3 = 1/XX has dimensions of inverse volume, 
on the volume of the system V and on the number of Since the term 'free volume' has other associations in the 
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Fig. 1. The geometric construction around vertex v in two 
dimensions. The vectors Vcv connect midpoints of ceU walls 
clockwise around each vertex. The vector Rev extends from the 
centre of triangle v to the centre of cell c. The antisymmetric 
part of each term in the tensor Cv ~ Rcvfcv describes 
the area of the quadrilateral of which Rev and rev are the 
diagonlas. This quadrilateral and a neighbouring one are shown 
shaded. The quadrilateral are exactly adjacent and therefore a 
sum over these terms give exactly the volume of the system. 

literature let us term here the analog of the free energy 
the free porosity, Y = — In Z/ (3. In direct analogy with 
thermodynamics, the (dimensionless) volumetric entropy 
is S" = P^dY/df] and the mean porosity is (V) = Y + S/p. 
Note that whether the ensemble is canonical or micro- 
canonical depends on whether the volume function is con- 
stant across the ensemble. The particular choice of ensem- 
ble is not essential to the following analysis and we shall 
not concern ourselves with it any further. 

The main difficulties in the application of this for- 
malism have been the identification of a volume function 



which is exactly additive when summed over the basic el- 
ements of the systems (cells or grains) and the isolation of 
a manageable phase space of independent degrees of free- 
dom. In the following we focus on the canonical ensemble 
at a fixed, but large, number of cells. The extension to 
grand canonical ensembles where cells are continually cre- 
ated and destroyed is straightforward. 

The main aim of this paper is to formulate a three- 
dimensional exact volume function in terms of local struc- 
tural descriptors and to demonstrate its use. However, as 
will become clear below, an important step toward this 
goal is the adaptation of a recent formulation of an exact 
volume function in two-dimensions 0]. 

2 Two dimensional structures 

A planar dry foam comprises of iV 3> 1 vertices, each of 
which connects three edges (trivalent structures). Vertices 
that connect more edges are rare in many common cir- 
cumstances and we shall disregard those in this paper. 
Nevertheless, the formalism developed here is quite gen- 
eral and applies to non-trivalent foams as well. The edges 
can be either straight or curved; the structural characteri- 
sation described here works for both cases and for systems 
that comprise of a mixture of such edges. The edges parti- 
tion the plane into N /2 polygonal (or polygonal-like in the 
case of curved edges) cells with an average of six edges per 
cell Uj and altogether 3iV/2 edges. This neglects boundary 
corrections which are at most of order 0{^/N). 

Consider the following network 1^. Around each ver- 
tex draw a triangle whose corners are the midpoints of 
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the edges that lead to it (see figure 1). We regard the the application of the compactivity-based formalism which 

sides of each triangle as vectors that circulate clockwise were mentioned in the introduction. The two-dimensional 

around the vertex. Consequently these vectors form closed volume function is 
loops around cells which circulate anticlockwise. A vec- 



(3) 



tor around vertex v and within cell c is indexed rev The yy^ — Ixr Cc 

L cv 

vectors rev form a network that spans the system, the r- 

We now wish to determine the phase space, i.e. the in- 

network. Let be the mean position (the centroid) of 

dependent degrees of freedom. The volumes Acv are ex- 

the corners of the triangle around v and Pc be the cen- 

pressed in terms of the vectors r^y and Rev and we note 

troid of the midpoints of the edges surrounding cell c. We 

that there are altogether 3iV vectors of each. But these 

define a vector that extends from p„ to pc, Rev ~ Pc~ Pv 

are not all independent. First, the self-duality between 

The network of R vectors is dual to the r-network. Each 

the R and r networks means that all the R vectors can 

rev-Rev pair forms a quadrilateral of which they are the 

be determined in terms of linear combinations of the r 

diagonals, such as the two shown shaded in figure 1. The 

vectors. Specifically, every Rev can be written as a lin- 

shape and geometry of a quadrilateral is characterised by 

ear combination of exactly {Ze + Zv — 2) r vectors, where 

the outer product 

Ze is the number of edges of cell c and Zv is the number 
of vectors around vertex v (in trivalent cellular systems 
™ '^^ Zy — 3, but it is arbitrary in granular assemblies and in 

where i,j = x,y are Cartesian components. The tensor non-trivalent foams). To see this recall that Rev extends 
Cv — J2c ^cvi where the sum runs over the three cells sur- from the centroid of vertex u, p^, = -J- Y^'^^^i {Zv — c')reiv, 
rounding vertex w, has been found to play a central role to that of cell c, pe — ^ X]f'=i^(^c ~ v')rev' ■ Therefore, 
in the analysis of stress transmission both in cellular sys- Rev ^ Pc Pv = ^ in^n can be expressed as a 

tems |S] and in granular assemblies 7 . The antisymmet- linear combination of the vectors circulating triangle v and 
ric part of Cev is Cev — C*^ = Aevi where e — (°jq) and cell c. An example is shown in figure 2 where i?3, which 

area of the quadtrilateral. The points to the centroid of a pentagonal cell, is a function of 
choice of the directions of r^v and Rev ensures that all the 5 + 3— 2 = 6r vectors, 
areas Aev are pseudo- vectors that point in the same direc- 
tion, into the page. This direction is taken to be positive. A ^ ^ 

Rs^T (4»^3 + 3r4 + 2r5 + rg) - - (2r2 + n) (4) 

key observation is that the quadrilaterals perfectly cover 

the plane, Asys = X^c-u^cu, providing an exact volume On average there are six edges around cells and therefore 
function. This volume function resolves the difficulties in an R vector depends, on average, on (Zv + Ze — 2) = 
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Fig. 2. The vector iia of the dual network is a linear com- 
bination of the contact network vectors r„ {n — 1 — 6). Note 
that it does not depend on the vector and therefore Rev 
depends only on (zc + 1) r vectors, pc is the position vector 
of the centroid of the cell (marked by X) . p„ is the centroid of 
the triangle made by the vectors ri, r2 and rr. 



7 r vectors. Thus, there remain only the 3N r vectors 
to examine for independence. These are also inter-related 
due of the loops that they form. The inter-dependence is 
directly evaluated in terms of the number of irreducible 
loops, namely, the elementary loops in terms of which all 
other loops can be expressed. In a loop of n vectors only 
n — 1 are independent, so every irreducible loop provides 
one dependent vector. There are only two types of such 
loops: the triangles around the vertices, of which there 
are N, and the polygons around the cells, of which there 
are N/2. This means that of the 3N r vectors, 3N/2 are 
dependent. The independent degrees of freedom are the 



components of the remaining 3N/2 vectors and therefore 
the phase space is 3iV-dimensional. 

We now note a surprising coincidence: the number of 
independent degrees of freedom is equal to the number of 
quadrilaterals! This, together with the fact that the sum 
over the areas of the quadrilaterals is the volume function, 
suggests that the quadrilaterals are the actual elementary 
'quasi-particles' of the system, not the vertices or the cells 
(nor the grains in the analogue context). It is convenient to 
express the phase space in terms of their areas, Acv = Aq. 
We term these elementary quasi-particles 'Quadrons'. In 
terms of the quadrons the partition function is 



(5) 



and the function 0{{Aq\) can be interpreted as a density 
of volumes (the analogue of the density of states in con- 
ventional statistical mecanics). It should be noted that 
cellular and granular systems support different limits of 
integration in In the former the lower limit on Aq 
may be as vanishingly small as the thickness of a cell wall 
and the upper limit, subject to the constraint that the to- 
tal number of vertices remains N , may be almost as large 
as the entire system. In granular packings Aq depends 
strongly on the grain size distribution and is bounded by 
finite upper and lower values due to the fact that grains 
cannot inter-penetrate. 

To illustrate the method, let us consider very large 
skeletal cellular systems (i.e. of vanishing cell wall thick- 
ness) and assume for simplicity that the quadrons are in- 
dependent, 0[{Aq}) — ■K^^{P{Aq). We are not aware of 
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any first-principles calculation of ©({Ag}) for disordered 
cellular structures and therefore we shall postulate several 
simple forms and analyse briefly the insight that they of- 
fer. First, let us consider the case of a uniform distribution 
between lower and upper volumes £ind Ai, repectively, 



(6) 



r ^ iiAo<Bg< Ai; 

{ otherwise , 
where 5 A = {Ai — Ao)/2. In dynamically evolving sys- 
tems 5 A and A = {Ai + Ao)/2 would scale with time. 
Substituting this form into the partition function gives 







3N 


z = 


/ e-P^dA 






Jao 





and 

S = 3N 



1 + ln 



sinh(/3M) 



(3e-l^'^smh{(3SA) 
]3SA 



(3SA 



3N 



(7) 



f36A tanh(/3(5A) 
The mean volume (or porosity) of this system per quadron 



is then 



{V)^A + l/f3 - 



5A 



ta.nh{pSA) 
and the mean porosity fluctuations is 



sinh^ipSA) 



As a second example consider a binary mixture of two 
types of independent quadrons, Aq and Ai, occuring at 
respective probabilities p and 1 — p 

ViA,) = p6iA, - Ao) + (1 - p)6{A, - A,) . (8) 
The one quadron partition function is found to be 



where u — ln[j3/(l — p)] and A and SA are defined as in the 
previous case. The entropy per quadron can be calculated 
from this expresion using s = S/N^of — Inz — pd\iiz/d(3 
and the mean porosity (volume) per quadron is 

{V) = A- Mtanh(/3M + u) . 

The mean of the porosity fluctuations per quadron is 

2 



SA 



^cosh{pSA + u), 
It is worth noting that as a function of the agitation the 
system can occupy volumes between A + SA and A + 
SAta.nh{u) and that this case is the exact analogue of the 
Bragg- Williams description of alloys [S]. 

For a third example we are prompted by observations, 
e.g. in soap froths TT, to consider a system with a broad 
distribution of cell sizes. Such distributions have generi- 
cally algebraic tails of the form 

ViAg)^CA'^-' ; 0<A,<a, (10) 

where < a; < 1 and C is a normalisation constant. The 
partition function is straightforward to compute 

Z=[Cp-^jix,Paf , (11) 

where 7(0;, y) is the incomplete gamma-function |^. As in 
the previous cases, the entropy can be obtained directly 
from this expression. The mean porosity per quadron is 

and the mean porosity fluctuations is 

^2, ^{x, I3a)^{x + 2, (3a) — ^f{x + 1, (3aY 



{SV 



2Vp(1 -p)e"'^^cosh {(35 A + u) 



(9) 



f3'^-i{x,l3af 
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Vertex v present the following discussion in the context of open-cell 

structures. 

Most common three-dimensional cellular structures are 
quadrivalent, namely, when each vertex connects to ex- 
actly four neighbour vertices. The skeleton of the struc- 
ture is obtained by extending lines between neighbouring 
vertices. The lines can be either straight or curved and 
represent the struts thinned down to vanishing widths. 
The skeletal struts meeting at a vertex form angles that 
are presumed in the following to be random. If they are 
Fig. 3. The geometric construction around vertex v. Vectors not (as may happen, e.g. due to dynamics processes dom- 
such as n and r2 connect midpoints of struts around each inated by surface tensions) then the configuration space 
vertex to form a tetrahedron. The vectors Rev points from jg spanned by fewer independent degrees of freedom. For 
the center of tetrahedron v to that of its neighbouring cell simplicity, wc shall not discuss this class of materials here, 
c. 5e„ = ri X r2/2 is the directed area of the face of the analysis can be extended to it straightforwardly, 

tetrahedron between v and c. 

Consider the following construction. Connect the mid- 
_ _, ,. . , points of the struts around every vertex by imaginary 

3 Three dimensional structures 

straight lines, which we term edges. The edges enclos- 
ing each vertex form a tetrahedron, as shown in figure 
We can use the insight from the two-dimensional analy- 3. Tetrahedra around two neighbouring vertices make con- 
sis as a basis for the formulation of the three-dimensional tact at the midpoint of the strut between them. The tetra- 
case. In three dimensions systems may be either close- hedra form an interconnected framework that spans the 
cell, in which case the cells are enclosed by material mem- system which we term, following intuition from granular 
branes (cell walls), or open-cell, in which case vertices are systems, the contact network. Every tetrahedron exposes 
connected by struts of given thicknesses and the cells are one triangular face to each of the four cells around it. De- 
interconnected. A sponge is an example of an open-cell noting a vertex by v and a cell by c, every such triangular 
structure. The formalism presented here is based on the face can be indexed uniquely by 'cu'. The triangular face 
structure of the skeleton (see below) and is therefore a de- that tetrahedron v exposes to cell c is characterised by 
scription of its topology. As such it applies equally well to Scv = Scvficv, where Scv is the area of the triangle and 



both types of structures. Nevertheless, it is convenient to ricv a unit normal to its plane, which points away from 
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the vertex and into cell c. The edges that make this trian- 4b. Rcvp and the edge of the triangle that it crosses define 

gular face arc regarded as vectors circulating the triangle a non-planar (skew) quadrilateral, an example of which 

in the clockwise direction when viewed from the vertex, is shown shaded in figure 4d. This and similar quadrilat- 

for example vectors ri and r2 in figure 3. erals tile the surface of the cell without overlap, leaving 

The centroid of the triangle Scv, Ccv, is defined as the ^o gaps, exactly as in two-dimensional systems. The only 

mean position of its three corners. We also define the cen- difference here is that the two parts of the quadrilateral, 

troid of tetrahedron v, py as the mean position of its four ^^e one within the triangular face and the one within the 

corners and the centroid of cell c, pe as the mean position Polygonal facet (shaded blue in figure 4) , lie in planes that 

of the midpoints of the struts that surround this cell. Note are tilted to one another due to the curvature of the sur- 

that while p„ is expected to be in the vicinity of vertex v ^^ce. Such a quadrilateral can be uniquely indexed by the 

it need not coincide with it. Nevertheless, for simplicity, cell, c, on whose surface it resides, by the polygonal facet, 

we shall refer to the point p„ from now on as vertex v ^'' surface and by the vertex, v, that gives rise to 

unless otherwise stated. original triangle Scv 

If we now extend straight lines from the four corners 



Let us inspect the contact network from inside cell 
c. Prom such a point the 'surface' of the cell consists of 
the triangular faces of the tertrahedra Scv which form an 
interconnected two-dimensional network. This structure is 
topologically identical (albeit on a closed surface) to the 
networks of triangles described above in planar systems. 
In particular, the triangles on this surface also enclose 
polygonal facets. 

The surface of the cell can now be tiled in exactly the 
same way that we did planar structures. We define the 
centroids of the non-triangular polygonal facets on the 
cell surface as the mean position of their corners and name 

^ cvp 

such a centroid (p. We extend three lines from the centroid Thus, we have arrived at the first goal - the parti- 

of every triangular facet, (cv, to the triangle corners (see tion of the volume of a three-dimensional cellular structure 

figure 4). Prom (cv to Cp we extend a vector Rcvp- This into volumes of elementary building blocks; the octahedra. 

vector crosses a particular edge, from the ends of which we These are the quadrons of the three-dimensional structure 

extend two lines to (p. This construction is shown in figure for the purpose of the entropic analysis Vcvp Vg. 



of such a quadrilateral to both the centroid of the cell, 
Pc and to the centroid of vertex pv, we get a non-convex 
octahedron (see figure 4f). The faces of this octahedron 
are shared with the faces of similar octahedra constructed 
on nearby quadrilaterals. The volume of the octahedron 
can be expressed in terms of the vectors Rcvp, Rev, and 
the appropriate edge of the triangle Scv, which we term 
rcvp- These octahedra cover the entire three-dimensional 
space without ovelaps, 
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Let US count the total number of quarons. A tetrahe- independent degrees of freedom. However, each facet is 

dron surrounding a vertex has four triangular faces, each common to two cells and therefore, when summing this 

giving rise to three quadrilaterals. Thus, with each vertex quantity over cells, every facet is exactly double-counted, 

we can associate a unique volume, composed of the vol- This gives that the total number of dependent vectors is 

umes of the twelve three-dimensional octagonal quadrons, "'//'^- Subtracting this number from the total number 

of struts, 2N, we obtain that the total number of inde- 

— Y^^^ _ (23) pendent degrees of freedom is 

cp 

The twelve quadrons around every vertex make a star-like 

Ndof^3{2N-Y,n}/2). (14) 
(stellated) icositetrahedron (i.e. a polyhedron of 24 faces) c 

J n J- J- • i 1 r ii -1 This sum can be evaluated using Euler's relation for the 

around it, with six laces extending into each oi the neigh- 

bouring cells of the vertex. We could then view the volume topology of the surface of a cell, 
of the system as covered by the N vertex icositetrahedra. 

Next, let us determine the dimensionality of the config- "-t + ~ ^ ' ^) 

urational phase space. As in two-dimensions, correlations together with the observations that in quadrivalent open- 
can only originate from vectors that close irreducible loops cell structures [TT] 
and our task is therefore to count these. The irreducible 

loops originate from all the polygonal facets on cell sur- — and ^^t — 6A^ . (16) 

c c 

faces, including the triangles. All other loops can be ex- ■ ,^ i • c n j_i ^ u r 

' ° ° Combining these relations gives nnaliy that the number ot 

pressed in terms of these. Consider then the surface of an , r r i 

^ degrees oi freedom is 

arbitrary cell c, consisting of vertices and facets. 
Because on the cell surface the vertices are trivalent (the 



fourth strut is part only of neighbouring cells) , the surface 



2A^-i5](2 + <-nS) 



3{N-N,) . (17) 



contains a total of = 3n^/2 struts. Each one of these Unlike in two-dimensions, it is not possible to simplify this 

struts can be described by a three-dimensional vector and result further to obtain Ndof in terms of A^ alone; Euler's 

so its specification requires three degrees of freedom. Ev- relation and the condition of quadrivalency are not suffi- 

ery vertex is quadrivalent and every strut ends in two ver- cient for this. Nevertheless, relation H17|) allows us to make 

tices, which gives that there are altogether 2A^ struts in two important observations. One is that Ndof < and 

the system. The n'j irreducible loops on the cell surface therefore Ndof is smaller than the number of quadrons. 

give n'j: dependent strut vectors on this surface, which in This means that, unlike in two dimensions, the quadrons 

turn gives that on this surface there are 3(3r7,J/2 — n'j) are inter-correlated, which has significant implications for 
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Fig. 4. The construction of a quadron: (a) Three tetrahedra on the boundary of one cell. (b,c) The two stages of the construction 
of a non-planar (skew) quadrilateral on the surface of the cell; First extend a triangle into the polygonal facet between triangles 
and then extend a triangle to the centre of the triangle Scv ■ The diagonal of the skew quadrilateral (red line in 4d and 4e) is a 
section of the line where the planes of the polygon facet and Scv intersect, (d) Connecting the four corners of the quadrilateral to 
the centroid of tetrahedron v. (e) Connecting the four corners of the quadrilateral to the centroid of the cell, (f) The quadron is 
a non-convex octahedron extending between the centroids of the vertex and of the cell. There are twelve such octahedra around 
vertex v whose joint volume makes a non-convex (stellated) icositetrahedron (i.e. of 24 faces). This is the volume associated with 
vertex v. These icositetrahedra perfectly fill the three-dimensional system and can also be used to partition its entire volume. 

the construction of the partition function. The second ob- three-dimensional cellular systems is even smaller than 
servation is that N^of < 3iV, indicating that contrary to that of two-dimensional systems, 
naive intuition, the dimensionality of the phase space of 



The partition function can now be written as 
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11 



Z = j0{{V,})Y[{e-^''^dV,), (18) Z^y^ Ce-^^'^-^^dV, 

(20) 

where 0{{Vq}) is the joint probabihty density of the vol- From this partition function we find that the entropy is 

umes Vg — Vcvp that are included explicitly in the integral. dlnZ 3(N — Nc)a? 

S ^ InZ ~ (3—— ^ (3"^ ^— ^ , (21) 

This probability density - the density of states - can be 

obtained from experimental measurements, from Simula- t^e mean volume (or porosity) per quadron is 
tions, or from analytic calculations on case-study systems. 

{V) = {Vo - Pa') , 

Armed with the form of this function it is possible to com- 
pute the entropy, the mean porosity, porosity fluctuations ^^^^ ^^^^ porosity fluctuations is 

and other related properties, in exactly the same way as , t, 

{SV') = [{{Vo - (ia'f + a')] • 

was outlined for two-dimensional systems following eq. 

For our second example we are inspired by data on 

Here too a first approximation would be to assume 

polymeric open-cell foams. In these materials the cell size 

that the Ndof quadrons are uncorrelated. This is the ana- 
distribution is generically skewed to the left We are 

logue of the ideal gas approximation and it simplifies the 

not aware of any fit to this distribution and therefore let 

computation of the partition function. To the best of our 

us assume that it can be described by the form 

knowledge there exists currently no data in the literature 
on the form of the density of states 0{{Vq]) in any sys- 
tem, either experimentally or analytically. Therefore, for ^(Kj) — CVq^^e~^''^^° ; <Vq < Vmax , (22) 
illustration purposes, let us first assume that the quadron 

with < a; < 1, < Vmax, Vq a typical size and C = 

volumes have a Gaussian density of states, 

Vq^^ /j{x, Vmax/Vo) a normalisation constant. For Vmax 

oo j{x, Vmax/Vo) — > r{x). It is plausible that the quadron 

^/x^x ^ (Vq-vpf volume distribution follows the same distribution as the 

V{Vg) = Ce . (19) 

cell sizes and, using the density of states (|^. the partition 
In this expression C is a normalisation constant and Vq function is 
fluctuates around a mean value Vq with a spread of a. If 
this is the steady-state distribution of dynamically evolv- 
ing systems then Vq and a would scale with time. Substi- 
tuting this form into the partition function we have 



ilix,Vmaxif3+l/Vo)) 



iVoP+ 1)^7 [X, Vmax /Va) 



Ndof 



(23) 
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Using now z — Z^l^'^°f we can readily compute the en- elementary entities are not vertices but rather quadrons, 

tropy per quadron s = Inz — d\nz/d(3, the mean porosity of which vertices are composed. In two dimensions and 

per quadron {V) — —d\nz/d[3, and the mean porosity in three dimensions a vertex is composed of three and 

fluctuations per quadron {5V'^) = 9^1nz/(9/3^. twelve quadrons, respectively, (iii) The volume functions 

The form given in H18|l may not be the most convenient make it possible to identify the basic variables on which 

form to represent the partition function. In some situa- the volumes depend, (iv) The functions make it possible 

tions it may be more convenient to use the more basic to determine the correlations between variables and so 

3(iV — Afc)-dimcnsional phase space of the components of identify the phase space of independent variables. We have 

the independent vectors. For example, when the aim is to found that in structure of N vertices and cells the 

compute expectation values of structural properties that phase spaces in two and three dimensions are 3iV and 

cannot be expressed directly in terms of quadron volumes. 3(A^ — N^), respectively, (v) The volume function in two 

Such properties could be the mean area of polygonal facets dimensions is derived directly from a fabric tensor which is 

on surface of cells, the mean moment of inertia of cells, or also instrumental for stress analysis in such structures, had 

the anisotropy of the cellular structure. In this case the they been in mechanical equilibrium. This may provide a 

volumes of the quadrons (or the vertices, for that matter) a bridge to a comprehensive theory of both stress and 

could be first expressed in terms of the vectors that con- entropy of cellular solids and solid foams. 

struct them, these expressions would be substituted into tj. ■ • ^ j.- ^ • ^ j- • ii j- 

' ^ it IS interesting that m two dimensions the dimension- 

the volume function W^, and finally a partition function t. r .t ^ ■ .^ ,i i r 

' J f ality oi the phase space is the same as the number oi 

of the form ITl could be written down and computed. , i -i • i- • • i n 

" quadrons, while m three dimensions it is much smaller. 

Moreover, we have found that the dimensionality of the 

4 Discussion phase space in three dimensions is smaller than that in two 

m ij 1 jij J.1 J.-J.- dimensions. The reason for this seemingly strange result 

io conclude, we have developed an exact characterisation ° ° 

r , 1 n i- 1 i r 11 1 i i ■ i is that, although the total number of variables in three- 

ot the configurationai entropy ot ceiluiar structures m two ° 

, ^, ,. . mi i 1- • 1 r 1 i- dimensional systems is much larger, so is the number of 

and three dimensions. Ihe two-dimensionai formulation ° 

, , J , ■ • • 1 i ■ i ii ii J- ■ 1 their inter-correlations. Therefore, the difference between 
has been used to gam insight into the three-dimensional 

mi 1 ■ ,1 , r i- -i J these two quantities, which gives the dimensionality of the 

case, ihe analysis uses the concept oi compactivity and i > t> j 

, i i 1 1 r i- 1-11 1 phase space, does not necessarily increase with the dimen- 

we have constructed volume lunctions, which have several i- ' 

, , /-N mi i • ii ii i 1- sion in which the physical system is embedded, 

advantages: (i) Ihey are exact m the sense that a linear ^ 

sum over elementary entities of the structures gives the Several case studies of densities of states, i.e. distribu- 

volume of the system, (ii) The analysis shows that the tions of quadron volumes, have been analysed. For these 
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cases we have calculated explicitly the mean porosity and 
the mean porosity fluctuations. The entropy has been cal- 
culated explicitly for one case in two dimensions and one 
case in three and its calculation has been outlined for the 
other cases. Since we have only aimed to illustrate the 
method, all the examples that we have studied contain no 
correlations between quadrons. It is not difficult to study 
such correlations for model systems - an example of such 
a calculation in two dimensions has been carried out in 
reference albeit for a granular system (but it applies 
equally well to cellular structures). Writing the partition 
function in terms of the density of states eliminates any 
explicit dependence on dimensionality, making the com- 
putation of the entropy in two- and three dimensions for- 
mally identical. Dimensional information is then contained 
implicitly both in the form of and in the number of 
quadrons. For ideal binary mixtures of quadrons we have 
derived an analogue of the Bragg- Williams description for 
alloys. The three-dimensional systems that we have pre- 
sented are idealised in that we have assumed no correla- 
tions between quadrons and therefore may be regarded as 
an 'ideal quadron gas approximation'. Note that the num- 
ber of independent quadrons is 3(A^ — Nc) - a small frac- 
tion (< 1/4) of the 12N quadrons in the system. Choosing 
these to be uniformaly distributed throughout the system 
means that their spatial density is quite low. The rest of 
the quadrons may be regarded as frozen degrees of free- 
dom. It is evident that the entropy depends sensitively on 
the density of states but, to our knowledge, there is 
no derivation in the literature of its form in any system. 



nor is there experimental or numerical data on it. How- 
ever, several works on cellular systems have analysed cell 
size distributions and by adopting these to quadron size 
distributions it may be possible to make progress. This is 
the route we have taken in our second example in three- 
dimensions. 
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